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Abstract 

We study the quantization of a classical system of interacting particles obeying a recently pro- 
posed kinetic interaction principle (KIP) [G. Kaniadakis, Physica A 296, 405 (2001)]. The KIP 
fixes the expression of the Fokker-Planck equation describing the kinetic evolution of the system 
and imposes the form of its entropy. In the framework of canonical quantization, we introduce 
a class of nonlinear Schrodinger equations (NSEs) with complex nonlinearities, describing, in the 
mean field approximation, a system of collectively interacting particles whose underlying kinet- 
ics is governed by the KIP. We derive the Ehrenfest relations and discuss the main constants of 
motion arising in this model. By means of a nonlinear gauge transformation of third kind it is 
shown that in the case of constant diffusion and linear drift the class of NSEs obeying the KIP is 
gauge-equivalent to another class of NSEs containing purely real nonlinearities depending only on 
the field p= |Vf. 
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I. INTRODUCTION 



A wide class of diffusive processes in nature, known as normal diffusion, are successfully 
described by the linear Fokker-Planck equation. Its relation to Boltzmann-Gibbs entropy 
p tr0py) in th e f _, 0f th e _ M e tnernte,^ is Well e, aMish ed flfl 

However, nonlinear Fokker-Planck equations f NFPEs) J, S, 0, Q, Q and their connection in 
the field of the generalized thermodynamics P, \v\ is nowadays an intense research area. 
In particular, many physical phenomena, in presence of memory effects, nonlocal effects, 
long-range effects or, more in general, nonlinear effects, are well understood with the help 
of NFPEs. 



To cite a few, let us recall the problem of diffusion in polymers {12], on liquid surfaces 
[lsI ]. in Levy flights [14] and enhanced diffusion in active intracellular transport Q|. Many 
anomalous diffusion systems have a quantum nature, like for instance charge transport in 
anomalous, joUd. Q, subrecoil laser coding Q and the aging effect in quantum dilative 
systems [18] . 

A still open question concerns the dynamics underlying the nonlinear kinetics governing 
the above anomalous systems. Langevin-like, Fokker-Planck-like or Boltzmann-like equa- 
tions have been used by different authors to generate nonlinear terms in the Schrodinger 



equation with the aim o 



particle interactions jlfll. 120 



describing, in the mean field approximation, the many quantum 
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It is now widely recognized that the presence of a nonlinear drift term as well as the presence 
of a diffusive term in a quantum particle current originates complex nonlinearities in the 
evolution equation for the ip-wave function. 

Different examples are known in literature of nonlinear Schrodinger equations (NSEs) orig- 
inating from the study of the kinetics governing the many-body quantum system. For 
instance, the Doebner-Goldin family equations [2^ have been introduced from topological 
considerations as the most general class of Schrodinger equations compatible with the lin- 



ear Fokker-Planck equation. In Ref. [24j the authors introduced a NSE starting from a 
generalized exclusion-inclusion principle (EIP) in order to describe systems of quantum par- 
ticles with different statistics interpolating with continuity between the Bose-Einstein and 



the Fermi-Dirac ones. In Ref. 25], in the stochastic quantization framework, starting from 
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the most general nonlinear kinetics containing a nonlinear drift term and compatible with 
a linear diffusion term, a class of NSEs with a complex nonlinearity was obtained. 

Recently, a kinetic interaction principle (KIP) has been proposed Q to define a special 
collective interaction among the iV-identical particles of a classical system. On the one 
hand, the KIP imposes the form of the generalized entropy associated with the system, 
while on the other hand it governs the evolution of the system toward equilibrium by fixing 
the expression of the nonlinear current of particles in the NFPE, thus governing the kinetics 
underlying the system. 

The link between the generalized entropic functional and the corresponding NFPE can also 
be obtained starting from a maximum entropic production principle. In Refs. ((J z| > taking 
into account a variational principle maximizing the dissipation rate of a generalized free 
energy, the authors obtained a NFPE in the Smoluchowski limit. The same NFPE was 
oUaled in M . fl fan a spastic process described by a generalized Langevin equation 
where the strength of the noise is assumed to depend on the density of the particle. 

In the present paper we perform the quantization of a classical system obeying KIP, where 
the statistical information is supplied by a very general entropy. 

Up to today, different methods have been proposed for the microscopic description of sys- 
tems. Schrodinger's wave mechanics, Heisenberg's matrix mechanics or Feynman's path- 
integral mechanics are some of the many. Another approach is given by the hydrodynamic 
theory of quantum mechanics originally owing to Madelung |2j| and de Broglie [28] and 
successively reconsidered by Bohm j^sf i n connection with his theory of hidden variables. 
In the hydrodynamic formulation of quantum mechanics, the complex linear Schrodinger 
equation is replaced by two real nonlinear differential equations for two independent fields: 
the probability density and its velocity field. Basically, such equations are formally similar 
to the equations of continuity and the Euler equation of ordinary hydrodynamics. 
This formalism is fruitful, as in the present situation, when the expression of the quantum 
continuity equation is inherited from the one describing the kinetics of the ancestor classical 
system. However, for a complete quantum mechanical description, besides the continuity 
equation, we need to know if and how we should generalize the Euler equation that describes 
the dynamics of the system. In this paper, in order to fix the nonlinear terms in the Euler 
equation, we require that the whole model be formulated in the canonical formalism. 
We obtain a class of NSEs with complex nonlinearity describing a quantum system of inter- 
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acting particles obeying the KIP in the mean field approximation. We study the case of a 
quantum system undergoing a constant diffusion process. The generalization to the case of 
a nonconstant diffusive process is also presented at the end of the paper. It is shown that 
the form of the entropy of the ancestor classical system fixes the nonlinearity appearing in 
;he evolution equation. By means of a recently proposed nonlinear gauge transformation 
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3l| this family of evolution equations is transformed into another one describing a 



nondiffusive process. In particular, when the kinetics of the system is governed by a linear 
drift term, the new family of NSEs contains a purely real nonlinearity depending only on 
the density of particles p = \ip\ 2 - 

As working examples we present the quantization of some classical systems described by 
entropies already known in the literature: BG-entropy, Tsallis-entropy j^], Kaniadakis- 
entropy j3] and the interpolating quantum statistics entropy 3^|. 

The plan of the paper is the following. In Section II we recall the relation between a 
given generalized entropy and the associated NFPE describing the kinetic evolution of the 
classical system in the nonequilibrium thermodynamic framework. This kinetic equation is 
justified on the ground of KIP. In Section III, firstly first present an overall summing up 
of the hydrodynamic formulation of the linear Schrodinger equation, then we generalize the 
method to quantize the classical system obeying EIP. The Hamiltonian formulation of this 
model is presented and a family of NSEs with complex nonlinearity is obtained. In Section 
IV we study the Ehrenfest relations and discuss the conserved mean quantities. In Section 
V, the nonlinear gauge transformation is introduced. Some relevant examples are presented 
in Section VI. The final Section VII present comments and conclusions. In Appendix A we 
give the derivation of the Ehrenfest relations while in Appendix B we briefly discuss the 
generalization of the model for a quantum system whose kinetics undergoes a nonconstant 
diffusive process. 



II. NONLINEAR FOKKER-PLANCK EQUATION 



Our starting point, according to nonlinear kinetics, is to relate the production of the 



entropy of a classical system to a Fokker-P 
following the classical approach to diffusion 



anck equation. This can be accomplished by 

J, El. 



We start by assuming a very general trace-form expression for the entropy (throughout this 
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paper, we use units with the Boltzmann constant k B set equal to unity) 

S(p) = - J dx J dp ln«(p) , (2.1) 

where n(p) is an arbitrary functional of the density particles field p = p(t, x), with x = 
(x 1 , ■ ■ • , x n ) a point in the n-dimensional space. 
The constraints 

Jpdx = l, (2.2) 

on the normalization and 

Js{x)pdx = E, (2.3) 

total energy of the system, with £(x) = p 2 /2m + V(x) the energy for each particle, are 
accounted for by introducing the constrained entropic functional 

S{p) = - J dx J dp \nn(p)-P J £(x)pdx- (3' J pdx . (2.4) 

The two constants (3 and /3' are the Lagrange multipliers associated with constraints (|2.2|) 
and Q . 

Quite generally, the evolution of the field p in the configuration space is governed by the 
continuity equation 

|^ + V-J = 0, (2.5) 

with V = (d/d x 1) • • • , d/dx n ), and assures the conservation of the constraint ()2.2j) in time. 
We assume a nonlinear relation between the current J and the constrained thermodynamic 
force 

Hp) = v g) , (2.6) 

by posing 

J = D 1 (p)F( P ) , (2.7) 

with D the diffusion coefficient and 7(p) still an arbitrary functional of p. 

Putting Eq. (j2.7J) in Eq. ()2.5|) . and taking into account the expression of S given in Eq. 

()2.4|) we obtain the following continuity equation 



dp 



+ V- {-D 1 (p)v[(3S(x)+(3' + \iik( P )]} = . (2.8) 



dt 

Introducing drift velocity 

u driit = -D(3V£(x) , (2.9) 



Eq. (JUBJ) takes the form of a NFPE for the field p 

+ V • [u drift7 (p) - Df(p) Vp] = , 



(2.10) 



where 



f(p) = 7(P) 



9 In k(p) 



(2.11) 



<9p 



Total current J = ./drift + Jdifr is the sum of a nonlinear drift current ./drift = w drift7(p)> 
and a nonlinear diffusion current Jdifr — ~Df(p) V p, different from Fick's standard one 
«/Fick = —D V p, which is recovered by posing j(p) = n{p) = p. 

Eq. (J2.1(J|) describes a class of nonlinear diffusive processes varying the functionals j(p) and 
K(p). 

We observe that for any given entropy (|2.1|) an infinity of associated NFPEs exists, one for 
any choice of Mo) . 

In Refs. jfllzj], starting from a variational principle which maximizes the dissipation rate 
of a generalized free energy functional, substantially equivalent to Eq. ()2.4jl . a NFPE in the 
position space as in Eq. (|2.1()jl has been obtained. The same NFPE fl2.1()j) was also obtained 
in Ref. |8( , starting from a stochastic process described by a generalized Langevin equation, 
where the strength of the noise is assumed to depend on the density of the particle. The 
nonlinear current, as in Eq. (|2.7|) . is given by the gradient of the functional derivative of a 
generalized free energy equivalent to Eq. (|2.4jh 

In Ref. jj] the problem of the NFPE derived from generalized linear nonequilibrium ther- 
modynamics was also discussed at length. 

At equilibrium, the particle current must vanish, and from Eq. ()2.6|) it follows 



where, without loss of generality, we posed the integration constant equal to zero (otherwise 
it can be included in the Lagrange multiply [3'). 
We obtain the equilibrium distribution of the system 



In particular, with the choice n(p) = e p, Eq. ()2.1)1 reduces to standard BG-entropy and 
Eq. (I2.13J) gives the well-known Gibbs-distribution. 



\nK( Peq )+pS(x)+P' = , 



(2.12) 




(2.13) 
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Let us now justify Eq. (j2.10j) starting from the kinetic approach introduced in |26J] through 
the KIP. In accordance with the arguments presented in Ref . [2^ , we consider the following 
classical Markovian process 

|f = / Ht, V^x)- 7r(t, x -> y)] dy , (2.14) 

describing the kinetics of a system of iV-identical interacting particles. 

For transition probability ir(t, x — > y) we assume a suitable expression in terms of the 

populations of the initial site x and the final site y. 

According to KIP we pose 

7r(t, a; -> y) = r(t, x, x - y) j(p, p) , (2.15) 

where p = p(t, x) and p' = pit, y) are the particle density functions in the starting site 
x and in the arrival site y respectively, whereas r(t, x, x — y) is the transition rate which 
depends only on the starting x and arrival y sites, during particle transition x — > y. 
The functional y(p, p') can be factorized in 

1 (p,p') = a(p)b{p')c(p,p') . (2.16) 

The first factor a(p) is a functional of the particle population p of the starting site and 
satisfies the boundary condition a(0) = 0, since if the starting site is empty transition 
probability is equal to zero. The second factor b(p') is a functional of the particle population 
p' at the arrival site, and satisfies the condition 6(0) = 1, because the transition probability 
does not depend on the arrival site if particles are absent there. Finally, the third factor 
c(p, p') takes into account that the populations of the two sites can eventually affect the 
transition collectively and symmetrically. 

The expression of the functional b(p') plays a very important role in the particle kinetics 
because it can stimulate or inhibit the transition x — > y, allowing, in this way, interactions 
originating from collective effects. 

With the assumptions made in Eqs. ()2.15|) and ()2.16j) for transition probability, according to 
the Kramers- Moyal expansion and assuming the first neighbor approximation, we can expand 
up to the second order the quantities r(t, y, y — x) j(p(t, y), p(t, x)) and j(p(t, x), p(t, y)) 
in Taylor series of y = x + u and y = x — u, respectively, in an interval around x, for u <C x. 
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We obtain 



r(i, x + u, u) j(p(t, x + u), p(t, x) 

\r(t, y, u)~{(p{t, y), p(t, xfj + [r{t, y, u)^p{t, y), p(t, x)) 



+ - 



1 d 2 



2 dy i dy j 



it, y, u)i(j)(t, y), p(t, x) 



uu. 

• J 



y^x 



and 



l(p(t, x), pit, x-uj) = I j(p(t, x), p[t, yfj - — j(p(t, x), pit, y)) u t 

i d 2 



+- 



2 dy.dy. 



l(p(t, x), pit, yfj u i u j 



y^x 



Using Eqs. (j2~TTj) and fll^ in Eq. (jZUSP , from Eq. (j2~TIjl it follows 



dp d 
dt dx. 



with i — 1, ■ ■ ■ , n and summation over repeated indices is assumed. 
InEq. fl2~P9l 

7(p) = 7(P, p') 



and 



a(p) 
Kp) 



while the coefficients and £ r are given by 



C = / r (*> w ) M * rfw ' 
^ = 2 / r((t ' V ' ^ Ul Uj dU ' 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



^2.22) 



^2.23) 



Defining (uj drift = — C ~~ the z-th component of thrift, and assuming the inde- 

pendence of motion in different directions of the isotropic configuration space we can pose 
= D6 ir It is easy to see that Eq. (I2~T9"|) reduces to Eq. (ETTUJ) . 

In conclusion we observe that Eq. (|2.1(Jj) is a NFPE in the Smoluchowski limit since it 
describes a kinetic process in the position space rather than in the phase space. This is a 
suitable form for the quantum treatment of the following sections. The passage from the 
NFPE in the phase space to the NFPE in the position space was rigorously elaborated in 
Ref. in the limit of strong friction, by means of a Chapman-Enskog-like expansion. 
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III. CANONICAL QUANTIZATION 



A. Quantization in the hydrodynamic representation 

In the hydrodynamic representation, the quantum mechanics formulation, can readily be 
obtained from the standard Schrodinger equation 

ifi §7 = -^ + y ^' < 3J) 

where V(x) is a real external potential. The complex field ip = ip(t, x) describing the 
quantum system is related to the hydrodynamic fields pit, x) and £(i, x) through polar 
decomposition j^, Q| 



if)(t, x) = p L/2 (t, x) exp f - E(t, x)j . (3.2) 
Eq. (|3.1|) is separated into a couple of real equations 

m || + m ( 9 .V)9 = v(£^-^0), (3.3) 
§? + V-i.-0, (3.4) 

where quantum velocity v, which in the linear case coincides with quantum drift velocity 
ttdrift, is related to the phase E(i, x) through 

m? = VS(i, x), (3.5) 

and 

3 = PV, (3.6) 

is the same relationship between current and velocity of the standard hydrodynamic theory. 
We remark that the quantum current ()3.6|) contains only a linear drift term. 
According to the orthodox interpretation of quantum mechanics the quantity p(t, x) = 
\ip{t, x)\ 2 represents the position probability density of the system normalized as 
/ p(t, x) dx = 1. 

Eqs. (I3.3j) - (|3.6|) form the basis of the hydrodynamic formulation which consists of a 
quasi classical approach to quantum mechanics. In this picture the evolution of the system 
can be interpreted in terms of a flowing fluid with density p(t, x) associated with a local 
velocity field v(t, x). The dynamics of such fluid is described by the Euler equation ()3.3|) 



and is governed by forces arising not only from the external field F ext {x) = — VV(a;) 
but also from an additional potential U q = — (h 2 /2m)A y /p/y i /p known as the quantum 
potential [2^. Remarkably, the expectation value for the quantum force vanishes at all 
times, i.e. (—V U q ) = 0. Finally, the continuity equation ()3.4|) assures the conservation of 
the normalization of wave function ip during the evolution of the system. 
Let us remark that the quantum fluid has a very special property. Because x) is a 
potential field for the quantum velocity, the quantum fluid is irrotational. As a consequence, 
in the linear Schrodinger theory, a non vanishing vorticity a;, defined by 

u) = V x v , (3.7) 

is possible only at the nodal region where neither £(£, x) nor VS(t, x) are well defined. 
At such a point V x V x) does not vanish in general, thus leading to the appearance 
of point- like vortices. 

Finally, putting Eq. ()3.5|) into Eq. (J3.3j) we obtain 

d£ (VE) 2 h 2 AJp . . 

— + - - ^=- + V(x = 0. 3.8 

dt 2m 2m v ; v ; 

This equation, in the classical limit h — > 0, reduces to the Hamilton- Jacobi equation for the 
function S. 

Eqs. ()3.4)1 and Eq. ()3.8j) can be obtained by means of the Hamiltonian equations 

as 5H 



dt 6p 
dp _ 5H 
~dt ~ III ' 



(3.9) 
(3.10) 



where the Hamiltonian 



with 



H = J H(p, T)dx , (3.11) 

JtoV = ^, + f<Z£+VV P . (3,2) 
2m 8m p 



represents the total energy of the quantum system. 



B. The many-body quantum system 



Let us now generalize the method described above by replacing the linear continuity 
equation Eq. ()3.4|) with the more general one obtained in analogy with the continuity 
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equation (|2.10|) describing the kinetics of a classical system obeying KIP. In the following we 
assume that the quantum system undergoes a constant diffusion process with D = const. 

We begin by introducing the wave function ip = ip(t, x) describing, in the mean field 
approximation, a system of quantum interacting particles. We postulate that the following 
NSE describes the evolution equation of the system 

^|y = -^A^ + A(r, ^ + V(x) V> , (3.13) 

where A(^*, if>) = W(ip*, if>) + % W(ip*, ip) is a complex nonlinearity, with W(ip*, ip) and 
W(ip*, ip) the real and the imaginary part, respectively. 

Using polar decomposition (|3.2|) . Eq. (|3.13|) is separated into a couple of real nonlinear 
evolution equations for phase and amplitude 

|f + + U q + W(p, E) + V(x) = , (3.14) 

^ + V-j o --pW(p, S)=0. (3.15) 

We require that both Eqs. ()3.14j) and ()3.15j) can be obtained through the Hamilton equations 
(I3.9|) - (j3.1()j) and, to accommodate nonlinearities W(p, S) and W(p, S), we introduce in the 
Hamiltonian density 7i an additional real nonlinear potential U(p, S) which describes the 
collective interaction between the particles belonging to the system 

n(p, E) = P + + U(p, E) + V(x) p . (3.16) 

2m 8m p 

By means of Eqs. (|3.9|) and (J3.1U|) it follows that the nonlinear functionals W(p, E) and 
W(p, S) are related to the nonlinear potential U (p, S) as 

W(p, H) = f U(p, E) dx , (3.17) 
o p 7 

W(p,E) = — — Ju(p,E)dx. (3.18) 

We assume that the quantum fluid satisfies a continuity equation formally equal to the 
classical one described by the NFPE (J27TUJ). By matching Eq. (I3~T5J) with Eq. (l2~TU|) we 
obtain the expression W and, accounting for Eq. (J3.18|) . we have the nonlinear potential 
U(p, E). Finally, the nonlinearity W(p, E), which follows from Eq. (|3.17|) . together with 
the quantum potential U q and the external potential V(x), describes the dynamic behavior 
of the quantum fluid according to Eq. ()3.14j) . 
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We observe that if the following equation holds 

J U(p, E) dx = V • F(p, E) , (3.19) 

with F(p, E) an arbitrary functional, taking into account Eq. (J3.18j) . Eq. (J3.15j) becomes 

d p 



_5_ 

5E 



dt 



+ V 



3 - F(p, E) 



. 



(3.20) 



Eq. ()3.19|) is fulfilled if functional U(p, E) depends on phase E only through its spatial 
derivatives j^ . 

Introducing the quantum drift velocity 



Adrift 



VE 



(3.21) 



which in the linear case coincides with the quantum velocity v given in Eq. ()3.5|) . and by 
comparing Eq. (J3.2(J|) with Eq. (J2.1(J|) we have 



F(p, S) 



VE 



m 



p- 7 (p) +Df(p) Vp. 



(3.22) 



By integrating Eq. ()3.18|) . the nonlinear potential assumes the expression 

(VE) 2 



U(p, E) 



2m 



7 (p)-p -Df(p)Vp-VZ + U(p) 



(3.23) 



where U(p) is an arbitrary real potential depending only on field p. Eqs. (13 .9)) and (J3.1(J|) 
give the following coupled nonlinear evolution equations 

as (VE) 2 dj(p) h 2 Ayfp 



at 2m op 2m ^fp \P ) 



(3.24) 



97 + v ' 



VE 



m 



1 (p)-Df(p)Vp 







(3.25) 



where G(p) = 5 J U(p) dx/5 p. 

In the classical limit h — > Eq. (|3.24j) becomes a nonlinear Hamilton- Jacobi equation for 
function E. It differs from the classical one owing to the presence of the nonlinear term which 
functionally depends on both p and E. We recall that such a nonlinearity was introduced 
consistently with the requirement of a final canonical formulation of the theory. 
We stress once again that in the approach described in this paper, we start from a nonlinear 
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generalization of the continuity equation that gives us only information on the kinetics. This 
equation is not enough to completely determine the time evolution of the quantum system. 
As a consequence, we have ample freedom in the definition of nonlinear potential U(p, E). 
Such freedom is reflected in the arbitrary functional U (p) which cannot be fixed only on 
the basis of the kinetic equation. There are many possible dynamic behaviors, one for any 
choice of U(p), compatible with the same kinetics. The nonlinear potential U(p) can be used 
to describe other possible interactions among the many particles of the system that have an 
origin different from the one introduced by the kinetic equation (|3.25|) . 

Actually, Eq. (|3.25jl is a quantum continuity equation for field p with a nonlinear quantum 
current given by 

3 = —-Y(p)-Df(p)Vp. (3.26) 
m 

We observe that, differently from the hydrodynamic formulation of the linear quantum 
mechanics, where the Bohm-Madelung fluid is irrotational, in nonlinear quantum theory the 
situation can be very different. In fact, by defining quantum velocity through Eq. (|3.6jh 
from Eq. (|3.2fi|) we have 

m $ = \E-mD InK(p)] , (3.27) 

p J 

which states the relationship between quantum velocity v and quantum drift velocity thrift 
for the nonlinear case. 

Expression ()3.27|) can be justified in terms of Clebsh potentials. In fact, as is well known, a 
nonvanishing vorticity can be accounted for in the Schrodinger theory by introducing three 
potentials p, v and A related to quantum velocity through the relation 

mv = V p + uV X . (3.28) 

Vorticity uj assumes a nonvanishing expression given by 

u = — Vz/xVA. (3.29) 
m 

By comparing Eq. ()3.28|) with Eq. (|3.27|) we readily obtain p = const, v = j(p)/p and 
A = £ — mD ln/c(p), respectively, and Eq. ()3.29|) becomes 

uf = —v( ^ ] xVE, (3.30) 
m \ P J 
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with no any contribution from the diffusive term. The irrotational case is recovered in linear 
drift 7(p) = p. 

The final expression of the NSE (|3.13|) is given by 

ih^ = -^Ai>+ [W{p, £) + i W(p, £)] V + ^(x) V , (3.31) 
with the nonlinearities 

W(p, E) = f M - l) g) 2 + m0/WV .(i) +GW , (3.32) 



and 



W(p, S) = - A V { [ 7 (p) - p] f 3 -A) + ^ V ■ [/(p) V p] . (3.33) 



2p r ,vr/ v p / J 2 p 



where a family of NSE was 



Eqs. (I3.32|) - (j3.33j) differ from the one obtained in Ref. 
derived in the stochastic quantization framework starting from the most general nonlinear 
classical kinetics compatible with constant diffusion coefficient D = h/2m. In particular, 
the real nonlinearity W arising in the stochastic quantization is found to depend only on 
field p, in contrast with expression ()3.32j) . where functional W depends on both fields p and 
E. 

Remarkably, we observe that when the kinetics of the system is governed by a linear drift, 
with j(p) = p, the expression of nonlinear terms ()3.32j) and ()3.33)1 simplify to 

W(p, E)=mD f(p) V • + G(p) , (3.34) 

and 

W(p, E) = V • \f(p) In «(p) V p] , (3.35) 

A p 

where /(p) = p(d/d p) ln/t(p). 

They are determined only through functional «(p) which also defines the entropy ()2.1|) of 
the ancestor classical system. 

IV. EHRENFEST RELATIONS AND CONSERVED QUANTITIES 

In this section we study the time evolution of the most important physical observables 
of the system described by the Hamiltonian density (|3.16J) with the nonlinear potential 
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()3.23|) : mass center, linear and angular momentum and total energy. The proofs are given 
in Appendix A. 

Let us recall that, given an Hermitian operator = 0^ associated with a physical observable, 
its time evolution is given by 

\ 0) ^(( s JLo^- r o s JL\ ix + m, (4.i) 



dV ' hJ \Sip Sip* J \ dt 

where the mean value (O) — J tp*Oipdx. The last term in Eq. (j4.1J) takes into account a 
possible explicit time dependence on the operator O. 
Observing that the NSE (j3.31|) can be written in 

ih^ = HiP, (4.2) 

where 

h 2 

H = — - — A + W(p, S) + i W(p, S) + V(x) , (4.3) 
2m 

Eq. (j4.1j) assumes the equivalent expression 



d 



1 /r, IdO* 



dt <0> = 5 <[R.H,0]) + ; <{toH,0}) + ( 7n .), (4.4) 



h\t ' Jf \dt 



where [■, ■] and {-, •} indicate the commutator and the anticommutator, respectively. 

By choosing O = x, from Eq. (j4.1j) we obtain the first Ehrenfest relation for the time 
evolution of the mass center of the system 

v mc = ^(x} = (^- u drift ^ . (4.5) 

We observe that only drift nonlinearity appears in this equation whereas the diffusion term 
makes no contribution. Eq. ()4.5j) states that, quite generally, v mc is not a motion con- 
stant. This fact implies that the quantum system is not Galilei invariant. The origin of the 
nonconservation of be found in the difference between quantity p mc = m v mc and 

the expectation value of the momentum operator p = (—ih'V) = J p'VJ^dx. These two 
quantities are equivalent only in the linear drift case. Differently from the former, the latter 
is in all cases conserved during the time evolution of the system, in absence of the external 
potential. This can be shown by means of the second Ehrenfest relation which follows from 
Eq. (I4~T1) by posing O = -ihV 

d 



lt (p) = {F ext (x)) . (4.6) 
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The time evolution of the expectation value of momentum is governed only by external 
potential V(x). On the average, the KIP introduce no effect on the dynamics of the system. 
This is a consequence of the invariance of nonlinearity W[p, E] + % W[p, S] under uniform 
space translation. 

In the same way, accounting for the invariance of nonlinearity for uniform rotations, the 
third Ehrenfest relation follows 

d 



lt (L) = {M ext (x)) , (4.7) 

where M ext (x) = x x _F ext (£c) is the momentum of the external force field. Eq. ()4.7|) 
is obtained from Eq. ()4.1|) after posing O = x x (— ifoV). Again, the nonlinear terms 
introduced by KIP as well as nonlinearity G(p) make no contribution, on the average, to 
angular momentum. 

Finally, the last relation concerns the total energy of the system given by the Hamiltonian 
E = H. By posing 

= ~A + -U{p, E)+V(x) , (4.8) 
2m p 

we have (O) = E and from Eq. (|4.1|) we obtain 

dE . . 

7t=°- (49) 

In conclusion, for a constant diffusion process we have shown that in absence of the ex- 
ternal potential the system admits three constants of motion: total linear momentum (p), 
total angular momentum (L) and total energy E. Such conserved quantities, according to 
the Noether theorem, follow as a consequence of the invariance of the system under uni- 
form space-time translation and uniform rotation. Moreover, the system is also invariant 
for global U(l) transformation which implies conservation of the normalization of field ip 
throughout the evolution of the system. 

In Appendix B we briefly discuss the case of a quantum system with a diffusion coeffi- 
cient D(t, x) that depends on time and position. This space-time dependence destroys the 
invariance of the system under uniform space-time translation and space rotation. As a 
consequence, all quantities (p), (L) and E are no longer conserved, even for a vanishing 
external potential. 

It should be remarked that the results discussed here, although very general in that they 
are independent of the form of nonlinearities W and W, are valid only for the class of the 
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canonical systems. In literature there are many noncanonical NSEs, obtained starting from 
certain physically motivated conditions, which are worthy of being taken into account. For 
these equations, the expression of H appearing on the right hand side of the Schrodinger 
equation cannot be obtained from Eqs. (|3.9j) and (J3.1(Jj) by means of a Hamiltonian function 
H = JHdx. 

Despite this, even for these noncanonical systems the time evolution of the mean values 
of the quantum operators associated with the observables can be derived through Eq. 
f!4.4|) . but what is important is that these operators can assume a different definition with 
respect to the one given in the canonical theory. For instance, in the canonical framework 
the energy is supplied by the Hamiltonian H of the system, whereas in a noncanonical 
theory it is identified with the operator ihd/dt = H. (We remark that in the canonical 
framework H and H are, in general, different quantities). Moreover, for a noncanonical 
theory, conservation of the energy and the momentum do not follow merely from the 
principle of invariance of the system under space-time translation. Their time evolution 
depends on the expression of the nonlinearities appearing in the Schrodinger equation. All 
of this clearly causes a profound difference in the resulting Ehrenfest relations. 
For instance, in Ref. 20( a noncanonical Schrodinger equation with complex nonlinearity 
was derived starting from a Fokker-Planck equation for density field p by assuming some 
physically justified separability conditions. The resulting evolution equation has the real and 
the imaginary nonlinearity given by W(p, £) = 7 (£ — (£)) and W(p, X) = (hD/2) Ap/p, 
respectively, where 7 is a constant related to diffusion coefficient D and such that D — > 
if 7 — > 0. It is easy to see that such nonlinearities cannot be obtained starting from a 
nonlinear potential U(p, E) through Eqs. (j3.17J) and (j3.18J) . The system described by 
this NSE turns out to be dumped and dissipative, even in presence of a constant diffusive 
process. In fact, it can be shown that, following Ref. |20j], from Eq. ()4.4|) it follows 
d (p)/dt = (F cxt ) — 7 (p) and dE/dt = d (E)/dt = — (7/m) (p 2 ), which is a very different 
situation with respect to the one discussed in the present paper, with the exception of the 
trivial case 7 = 0. 
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V. GAUGE EQUIVALENCE 

We introduce a nonlinear gauge transformation of the third kind 



ip — ► <p = ip exp (--mD In n(p) 



(5.1) 



which, being a unitary transformation, does not change the amplitude of wave function 
1^1 2 = 1 01 2 = Pi an d transforms the phase S of the old field ip, into phase a of the new field 
according to the equation 

(j = E-mD 1hk(p) . (5.2) 
Consequently, the nonlinear current (|3.2fij) takes the expression 



3 — 3 = -zr i{p) ■ 



m 



(5.3) 



with only a nonlinear drift term. 

Let us observe that, at the classical level, the similar transformation 

tt^ift = u dri ft - D V In k(p) , 



(5.4) 



changes total current J — ► J' = u^ rift 7(p) into another one consisting only of a nonlinear 
drift term. 

Performing the transformation (|5.1J1 . Eq. (J3.31j) becomes 



at 2m L J 

where the new nonlinearities W(p, a) and W(p, a) are given by 



(5.5) 



7 pgH (7) 



A(p)Ap + / 2 (p) (Vp) 1 



+ C(p) 



(5.6) 



with jf = pVcr/m, 



AO) = 7(p) 



In /c(p) 

op 



h{p) = ■= 



and 



W(p, E) 



2 <9p 



^V{ [7 (p)-p](^ 



(5.7) 
(5.8) 

(5.9) 
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Eq. ()5.5|) is still a NSE with a complex nonlinearity due to the presence of the nonlinear 
drift term in the quantum current expression ([5.3)1 . 

Basically, both equations (|3.31|) and (|5.5|) are different NSEs describing the same physical 
system. This is a consequence of the unitary structure of the transformation (|5.1|) which 
implies that the probability position density for field if) and field assumes the same value 
at any instant of time j^J. 

In the case of j(p) = p expressions f!5.6|) and ()5.9|) can be simplified and the NSE ()5.5|) 
assumes the form 



ih— - = A0 + mD 

at 2m 



fi{p)Ap + h(p) (Vp) s 



4> + G{p) <p + V(x) <f) 



with 



h(p) = P 

Up) 



dp 

1 df x {p) 

2 dp 



In 



(5.10) 

(5.11) 
(5.12) 



which contains a purely real nonlinearity depending only on field p. 

We observe that although Eq. (|5.1|) transforms the nonlinear current into another one 
without the diffusive term, NSEs (|5.5|) and (|5.1U|) contain a dependence from on diffusion 
coefficient D. 

The NSE (|5.5jl is still canonical. It can be obtained from the following Hamiltonian density 

n 2 (y P y 



H(p, a) 
with nonlinear potential 

U(p a) - 



< Va) 2 



2 m 



P + 



• m p 



U(p, a) + V(x)p 



2 m 



lip) ~ P 



^fi(p) (Vpf + U(p). 



(5.13) 



(5.14) 



In this sense Eq. ()5.1j) defines a canonical transformation. 

In conclusion, let us make the following observation. Eq. ()5.5|) admits the following 
continuity equation 



9 P „ 



Vff 



m 



lip) 



. 



(5.15) 



A natural question is: what kind of NSE is obtained if we quantize a classical system obeying 
the continuity equation d p/dt + V ■ J' = with the method described above? 
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We easily have 



dcj) h 

% Ti -7— = — 

at 2m 



dp 



h 



Jo 

VP, 



1 * 



<j> + G(p) + V(x) cf) 



(5.16) 



where now p and a are independent fields representing the amplitude and phase of wave 
function (p. Eq. ()5.16|) can be derived through the Hamiltonian density ()5.13|) with nonlinear 
potential 



2 m 



7 (p)-p +U(p) . 



Potentials (|5.14jl and (|5.17jl differ for the quantity 



U(p) = U(p, a)-U x {p, a) 



m D 2 



MP) (Vp) ; 



(5.17) 



(5.18) 



which depends only on field p. This nonlinear potential U(p) does not affect the continuity 
equation and thus cannot be obtained starting directly from Eq. (|5.15|) . 



VI. SOME EXAMPLES 



To illustrate the relevance and applicability of the theory described in the previous sec- 
tions, we derive and discuss some different NSEs obtained starting from kinetic equations 
known in literature. In the following Section, for simplicity's sake we omit the arbitrary non- 
linear potential U(p) and focus our attention only on the effect yield through the potential 
introduced by the KIP. 



A. Boltzmann-Gibbs-entropy 

It is well known that when the many body system is governed by short-range interactions, 
or when interaction energy is neglecting with respect to the total energy of the system, the 
suitable entropic functional is given by the BG-entropy 

S BG (p) = - J p\n(p) dx . (6.1) 

This entropy arises from Eq. (|2.1j) by posing n(p) = ep with a(p) = ep and b(p) = 1. It is 
readily seen that j(p) = epc(p). 
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Among the many NFPEs compatible with entropy (j6.1|) we consider the simplest case of 
linear drift by posing c(p) = 1/e. Then the continuity equation ()3.25|) becomes the standard 
linear Fokker-Planck equation 

dp 



-£ + V-(j o -DVp)=0, (6.2) 



at 

whereas the evolution equation for the quantum system is given by the following NSE 

8 tb h 2 ( i \ h Ao 

ih^-= Af + mDV- P i) + i-D — ip + V(x)ip , (6.3) 

at 2m \P J 2p 

which is recognized as the canonical sub-family of the class of Doebner-Goldin equations 

parameterized by diffusion coefficient D. We recall that Eq. (j6.2j) was obtained in the 

quantum mechanics theory startin g fr om the study of the physical interpretation of a certain 

family of diffeomorphismin group [23]. 

By performing gauge transformation (|5.1)l . Eq. ()6.3j) becomes 



ih—- = Acb + mD 

at 2m 



Ap 1 fVp\ 2 



<f> + V(x)<f>, (6.4) 



p 2 V p , 

which was studied previously in [3^|. In particular, Eq. ()6.4j) is equivalent to the following 
linear Schrodinger equation 

i*|7 = -^-A X + V(*)x, (6.5) 
at 2m 



with % = h y 1 — (2 m D/h) 2 and field x is related to hydrodynamic fields p and a through 
the relation \ = p 1 ^ 2 exp(i<r/£). 

This appear to be an interesting result. By quantizing a classical system described by 
MB-entropy the standard linear Schrodinger equation was obtained. In this equation the 
nonlinear terms describing the interaction between the many particles of the quantum system 
are absent. This is in accordance with the general statement that MB-entropy is suitable 
for describing systems with no (or negligible) interaction among the particles. 



B. Generalized entropies 

In presence of long-range interactions or memory effects persistent in time, it has been 
argued that MB-entropy may not be appropriate in describing such systems. For this reason, 
many different versions of Eq. (j6.1|) have been proposed in literature. 
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Very recently, Ref. [36 
tion 



introduced a bi-parametric deformation of the logarithmic func- 



ln M ( x ) = ^ > (6-6) 



which reduces, in the (k, r) — > (0, 0) limit, to the standard logarithm: ln {00} (x) = lnx. By 
replacing the logarithmic function in Eq. (j6.1j) with its generalized version (|6.6|) . we obtain 
a bi-parametric family of generalized entropies 

5 {K ,,.}(P) = ~/p ln {K ,, } (p) dx » ( 6 - 7 ) 



introduced, for the first time, in Refs. |38l. I39 | . Remarkably, this family of entropies includes, 
as special cases, some generalized entropies, well known in literature, used in the study of 
systems exhibiting distribution with asymptotic power law behavior. Among them we can 
cite Tsallis-entropy |32j which follows by posing r = ±|k| 

S q (p)=J^-dx, (6.8) 



with q — 1 ± 2 \k\ and Kaniadakis-entropy [26(, for r = 

s w (p) = -J 2k dx - ( 6 - 9 ) 

Both these entropies, as well as other one-parameter deformed entropies, originated from 
Eq. ()6.7|) (3?]], can be employed to describe generalized statistical systems such as, for in- 
stance, charge particles in electric and magnetic fields j^, 2d-turbulence in pure-electron 



plasma 



43, 



44|. 



4l| . Bremsstrahlung 



42J and anomalous diffusion of the correlated and Levy type 



In addition to the many applications where Tsallis-entropy has been employed [45], 
Kaniadakis-entropy (|6.9|) has been successfully applied in the description of the energy dis- 
tribution of fluxes of cosmic rays [26], whereas the entropy in ()6.7|) with k 2 = (r + l) 2 — 1 
has been applied in the generalized statistical mechanical study of g-deformed oscillators in 



the frame-work of quantum-groups 



46] 



Despite the topics recalled above, there is currently great interest in studying quantum sys- 
tems with long-range microscopic interactions. Systems such as quantum wires, which are 
now possible in practice thanks to recent technological advances, require on the theoretical 
ground, the development of a quantum (nonlinear) theory capable of capturing the emergent 
facts |47J. 
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The entropy in (|6.7|) arises from Eq. ()2.1|) by posing 



ln/c(p) = A ln { - 



a 



(6.10) 



with A = (1 + r - K)( r+K )/ 2K /(1 + r + K )( r - K )/ 2K and a = [(1 + r - «)/(l + r + ft;)] 1 / 2 *. 
Among the many different possibilities, we discuss the case of linear drift with 7(p) = p. By 
taking into account Eq. (jfi.lOJ) we have continuity equation ()3.25j) with 



f(p) 



a + p r+K - a_p r K 



(6.11) 



where a± = (r±«)(l + r±/c)/2« are constants. 
The associated NSE assumes the expression 



ft 2 



z ft — — = — 
at 2m 



A</> + mL>' 



P 



f{p)Ap + f(p) (Vp) 2 + y(a;)0, (6.12) 



with 



f(p) = b + p r+K - l -b„p r -«- 1 , (6.13) 

and 5± = a± (r ± k — 1/2). 

Eq. ()6.12|) contains only a purely real nonlinearity and reduces to Eq. ()6.4|) in the («, r) — > 
(0, 0) limit, as well as Eq. (|6.7jl . which reduces to the standard BG-entropy. 
In particular, for Tsallis-entropy, the continuity equation ()3.25p . with 



f(p) =QP 



9-1 



(6.14) 



becomes the diffusive NFPE j4S| while the corresponding NSE is given through Eq. ()6.12j) 
with 



f(p) 



<i-- 2 )P 



ft- 2 



(6.15) 



and reduces to Eq. ()6.4|) in the q — > 1 limit just as entropy (j6.8|) reduces to BG-entropy. 
We observe that in Refs. j3, 3 the quantization of a classical system described by Tsallis- 
entropy has been already discussed. There, a NLS compatible with the continuity equation 
dp^/dt + V • (p A 'w drift ) = was obtained with a different approach. The nonlinearity 
appearing in the NLS of Refs. (3, |s(| reduces, for p = 1 and g — > 2 — g, to the same one 
reported here. 

On the other hand, for Kaniadakis-entropy, the continuity equation is given in Eq. ()3.25|) 
with 



f(p) 



1 



(6.16) 
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which coincides with that proposed in Ref. 
(l6~T2l with 



while the associated NSE is given in Eq. 



fa = h 



(K + D p k + (^-i) (« • \ ) 



(6.17) 

and reduces to Eq. ()6.4j) in the k — > limit just as entropy ()6.9|) reduces to BG-entropy. 



C. Interpolating bosons- fermions- entropy 

In Ref. 



33], on the basis of the generalized exclusion- inclusion principle the authors 
introduced a family of NFPEs describing the evolution of a classical system of particles whose 
statistical behavior interpolates between bosonic and fermionic particles. The equilibrium 
distribution governed by the EIP can be obtained by maximizing the following entropy 



Seip(p) 



p hip — — (1 + K p) ln(l + k p) 

K 



dx , (6.18) 



with — 1 < k < 1. In particular, for k = ±1 we recognize the well-known Bose-Einstein 
and Fermi-Dirac entropies, whereas intermediary behavior follows for — 1 < k < 1. Entropy 
(I6.18|) can be obtained from Eq. (|2.1|) by posing a(p) = p and b(p) — 1 + k p. 

Some examples of real physical systems where EIP can be usefully applied are to be found 
in the Bose-Einsten condensation. Typically, the cubic NSE is used to describe the behavior 
of the condensate by simulating in this way the statistical attraction between the many 
bodies constituting the system. In spite of the simplest cubic interaction, other interactions 
like the one introduced by the EIP can be adopted to simulate an attraction among the 
particles. 

In the opposite direction, almost-fermionic systems can be found in the study of the motion 
of electrons and holes in a semiconductor. In fact, while if separately considered electrons 
and holes are fermions, together they constitute an excited state behaving differently from 
a fermion or a boson. The same argument can be applied to the Cooper-pair in the super- 
conductivity theory. Such excitation differs from a pure boson state because of the spatial 
derealization of the two electrons, which are not completely overlying. Deviation from Bose 
statistics must be taken into account. 

In the following we discuss separately two different choices for functional 7(p). 
In the linear drift case, with c(p) = 1/(1 + up), the evolution equation for field p assumes 
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the expression 



dp 
dt 



D 



Vp 



which was proposed in Ref. 



k p 



(6.19) 



221 ] . By means of Eq. (|5.1|) . nonlinear current j 



D V p/(l + up) —>■ j assumes the standard bilinear form and the corresponding NSE fol- 
lows from Eq. (|5.1()|) 

2" 



i h 



30 
~dt 



2 m 



A0 + 



[1 + kpY 



Ap l-3np (Vp 
~~ 2(1 + k P ) \ p~ j 



<f) + V(x)(j) . (6.20) 



We can observe that in Eq. ()6.2())1 the EIP is accounting through a diffusion process and 
its effect vanishes in the D — > limit where it reduces to the standard linear Schrodinger 
equation. Eq. ()6.2U|) has a purely real nonlinearity depending only on field p. 

In a different way, by making the choice c(p) = 1, the continuity equation (|3.25|) becomes 



^ + V-[j o (l + K P )-DVp 



. 



(6.21) 



The gauge transformation changes nonlinear current j = j (1 + k p) — D V p — > j = 
j (1 + k p) containing only a nonlinear drift term and Eq. (|6.21j) reduces to 



dp ^ 



(6.22) 



This equation was introduced at the classical level in Ref. 33] and subsequently reconsidered 



at the quantum level in Ref. 



i n —— = — 

at 2m 



24j . The NSE associated with Eq. ()6.22j) is given by 

2" 



A6 + 



m 



k p 

2 



Ap 1 + 2/cp (Vp 



P 



2(1+ K P ) 



p 



+ K— ( - — 

p K^ + Kp, 



n 



i k — V • 

2p 



3 P 



<P + V{x)<f) . 



(6.23) 



1 + Kp / 

We observe that Eq. ()6.23|) still has a complex nonlinearity due to the nonlinear structure 
of quantum current j and both the nonlinearities W and W depend on fields p and a. 
Moreover, in Eq. ()6.23|) . EIP is accounted through a nonlinear drift term and survives even 
in absence of a diffusion process (D — > 0). 

Factor (1 + np) in nonlinear current j takes into account the EIP in the many particle 
system. In fact, transition probability ()2.15j) from site x to y is defined as n(t, x —>■ y) = 
r(t, x, x — > y) p(t, x) [1 + k p(t, y)]. For k ^ the EIP holds and parameter k quantifies 
to what extent particle kinetics is affected by the particle population of the arrival site. 
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If K > the 7r(£, x ^ y) contains an inclusion principle. In fact, the population density 
at arrival point y stimulates the particle transition and therefore transition probability 
increases linearly with p(t, y). Where k < the ir(t, x — > y) takes into account the Pauli 
exclusion principle. If the arrival point y is empty p(t, y) = 0, the ir(t, x — > y) depends 
only on the population of the starting point. If arrival site is populated < p(t, y) < p ma x, 
the transition is inhibited. The range of values that parameter k can assume is limited by 
the condition that 7r(t, x — ► y) be real and positive as r(t, x, x — > y). We may conclude 

that K > -1/pmax- 

A physical meaning of parameter k can be supplied by the following considerations. We 
recall that Bose-Einstein and Fermi-Dirac statistics originate from the fundamental principle 
of indistinguishability in quantum mechanics which is closely related to the symmetrization 
of the wave function. Completely symmetric wave functions are used to describe bosons while 
fermions are described by completely anti-symmetric wave functions. Thus, intermediate 
statistics arise in presence of incomplete symmetrization or anti-symmetrization of the wave 
function and the concept of degree of symmetrization or degree of anti-symmetrization has 
been introduced j^. Parameter k has the meaning of degree of indistinguishability of 
fermions or bosons, corresponding to the degree of symmetrization or anti-symmetrization, 
respectively. Value k — — 1 corresponds to the case of fermions and denotes a complete 
anti-symmetric wave function whereas value k — 1 corresponds to the case of bosons and 
denotes a complete symmetric wave function. In addition, value k = is associated with 
classical MB statistics and all the intermediate cases arise when k assumes all the values 
between —1 and 1. 

Eq. (|6.23p . for D = 0, was obtained previously in Ref. [24j, where the canonical quanti- 
zation of the classical system obeying EIP was accounted for. As discussed in Section El Eq. 
()6.23|) differs from the NSE obtained in 2M f° r a rea ^ nonlinearity originated from nonlinear 
potential U(p) = —mD 2 (V p) 2 / p (1 + Kp) and depending only on field p. 
Finally, we observe that different from Eq. (|fi.2()jl . Eq. (|fi.2H|) has vorticity different from 
zero. The Clebsh potentials corresponding to current j = (Vcr/m) p(l + up) are given by 
u — 1 + K p, X = a and p = const and vorticity assumes the expression 

w = -V«xV(T. (6.24) 
m 

In Ref. 0> localized, static, fermion-like vortex solutions (« < 0) were obtained and 
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studied starting from Eq. ()6.23|) with D = 0. We observe that in |5l|, |52j a different defini- 
tion of the Clebsh potentials corresponding to p = A = o and v = k p was adopted. Despite 
this, vorticity assumes the same expression that is given by Eq. (|6.24|) in both cases. 
EIP vortex solutions are important on the theoretical ground and for interpretation of ex- 
perimental results of several applications. For instance, they can be employed in the study 
of fermion-like vortices observed in 3 He-A superfluidity or in heavy fermion superconductors 
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53, 
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VII. CONCLUSIONS 

We have presented the quantization of a classical system of interacting particles obeying 
a kinetic interaction principle. The KIP both fixes the expression of the Fokker-Planck 
equation describing the kinetic evolution of the system and imposes the form of its entropy. 
In the framework of canonical quantization, we have introduced a class of NSEs with com- 
plex nonlinearity obtained from the classical system obeying KIP. The form of nonlinearity 
A(-0*, ip) is determined by functional n(p), which also fixes the form of the entropy of an- 
cestor classical system. 

Among the many interesting solutions of the family of NSEs (|3.31|) we observe that for a 
free system with V(x) = 0, and posing G(p) = 0, the planar wave 

ip(t, x) =Aexp(-^(Lut-k-x) \ , (7.1) 

with constant amplitude A = const is the simplest solution, where the relationship between 
uj and k is given by 



2 j_2 

(7-2) 

p=A 2 



nrhr dM P ) 

to 



2m dp 

and reduces to the standard dispersion relation for •y(p) = p. 

When the quantum system is in a stationary state such that d p s /dt = 0, the relationships 
between distribution p s and phase S s follow from Eq. ()3.25j) 

'•^-'H^-"'))' (73) 

which mimics the classical equilibrium distribution (j2.12)l . as can be seen by replacing 
T ls (x)/mD with —f3£{x). Despite this, we stress that such an analogy is purely for- 
mal. The equivalence between Eqs. f!2.12|) and ()7.3|) requires that the following relation 
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Y, s {x)/mD = —(3£{x) must hold. In the general case the expression of stationary phase 
£ s (a?) must be obtained from Eq. (|3.24jh after posing dT, s /dt = 0, with p given through 
Eq. (J73J). 

Finally, another interesting class of possible solutions are solitons. It is well know that soli- 
ton solutions in NSE arise when the dispersive effects, principally due to term —(h 2 /2m) A ip, 
is exactly balanced by the nonlinear terms. The existence of this class of solutions depends 
on the particular form of functionals 7(p) and n(p) which fix the expression of nonlinearities 
W(p, S) and W(p, E). A special situation, where soliton solutions are found within the 
NSEs derived in this paper, is given by the ElP-equation 1)6.23)1 with D = j^J where 
l{p) — P (1 + K P) an d k{p) = p/(l + ftp). 

The study of soliton solutions for other functional choices of 7(p) and n(p), like, for instance, 
the ones related to the generalized entropies discussed in Section VI-B, is a very important 
task which deserves further research. These solutions may lead to practical applications. In 
fact, in recent years there has been great interest in the formulation of models where solitons 
can interact with a long-range force Typical nonlinear models supporting solitons, like 
the sine-Gordon model, arise from short-range forces. However, there is experimental evi- 
dence that most real transfer mechanisms have long-range interaction, as noted in condensed 



matter theory 



571 ] or in spin glasses 
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APPENDIX A. 



We present proof of the Ehrenfest equations discussed in Section IIVI In the following we 
assume uniform boundary conditions on the fields in order to neglect the surface terms. 
Let us rewrite Eq. (j4.1j) in a more suitable form. Accounting for the relation 



_6_ 



'5 ih 5 
Jp~2p5E 



(A.l) 



Eq. (|4.1|) becomes 
d_ 



(O) 



5H 

Tp~ 



1 SH 



p 



O 



~dt 



(A.2) 



Eq. (|4.5|) can be obtained starting from Eq. ()A.2|) by posing O = x 
d 



dt 



(x) 



1p — Xlp — lp X— Tp 



Sp 



5p 



dx H — / 
2 J 
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p 5£ 



p 5E 



dx 
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x jk dx 

xV 
VE 



VE 



m 



1 (p)-Df(p)Vp 



dx 



rn 



1 (p)-Df(p)Vp 



dx 



I 7(P) dx-D / V F(p) dx 

J m J 

hip) 



' Wdrift 



(A.3) 



where 

F(P)= ff(p)dp'. 
To show the validity of Eq. (J4.6)) we pose O = — i h V in Eq. (jA.2|) so that 



5p 



5p 



da; 



dx 



/ ^ (V»*V^ + VVV*)-i^ / -^ (f Vi^-^Vf) dz 
J o p 2 J p oh 
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5 H 

V p + -rr^ V £ | dec 



where an integration by parts has been performed, and we have posed 



(A.4) 



(A.5) 



n 



Taking into account the relation 



VH= 5 -^-Vp+ 5 -^Vh + pVV(x) , 

dp o S 



from Eq. (|A.5|) it follows 



d 

It 



(p) = JvHdx- J pW(x)dx 



(F cxt {x)) 



(A.6) 
(A.7) 



(A.8) 



(A.9) 



Eq. ()4.7|) can easily be obtained following the same steps used in the proof of Eq. ()4.6|) . 
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Finally, by posing O = -(h 2 /2 m) A + U(p, S)/p + V(x) in Eq. (TPJ) . where C/(p, S) is 
given in Eq. ()3.23j) . we have 

5H ( h 2 . U \ , 17 ft 2 . C/ 
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dU Udp 
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dt p dt 

where a double integration by parts has been performed. Taking into account 

i h 



(A.10) 



2 m 

'ip* Aip + ip Atp* = 2 p 
which follow from Eq. (|3.2jl . and the relation 
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dU 
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Eq. (jA.lOJ) becomes 
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By using the relations 



1L 

2m 



VP 



h 
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T P 



5p 



m 



(A.15) 
(A.16) 



which follow from Eqs. (|3.11|) . ()3.16|) and ()3.23|) . and motion equations ()3.9|) and ()3.10|) . we 

obtain 
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APPENDIX B. 



We briefly discuss the generalization of the theory for quantum systems obeying the KIP 
and undergoing a diffusive process with a diffusion coefficient D(t, x) depending both on 
time and space position. 
Given the following Hamiltonian density 



H(p, S) 



fvs) 2 , , h 2 (v P ) 2 



2 m 



l(p) 



>m p 



D(t, x) 7 (p) V In K (p) ■ VS + U{p) + V{x) p 



(B.l) 



from the Hamilton equations (|3.9j) - (|3.1()j) we obtain the NSE 



% h 



dip 



h 7 



W(p, E) +iW(p, S) V + G(p) ip + V(x) * 



(B.2) 



dt 2m 
with nonlinearities 

^ ,S) = l(^7 _1 ) (j) + m 7(p)^ln^(p)V-^(t )a; )^+G(p),(B.3) 
and 

W(p, S) = -JL V {[ 7 (p) -/)]VS} + iv- [£>(f, aj) 7 (p) V ln*(p)] . (B.4) 
2m p 1 J 2p 

The system described by Hamiltonian (jB.l|) is dissipative since dE/dt ^ 0. This is a 
consequence of the time dependence of D which breaks the invariance of Eq. (jB.lj) under 
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uniform time translation. In the same way, linear momentum as well as angular momentum 
are no longer conserved, even in absence of the external potential, as a consequence of the 
position dependence of D which breaks the invariance of Eq. (|B.1J) under uniform space 
translation and uniform space rotation. This can also be seen from the Ehrenfest relations 



where A(p, S) = f(p) V lnp • w drift . 

Finally, the gauge transformation described in Section V cannot be performed, in general, 
when the diffusion coefficient has spatial dependence. In fact, the transformation in (|5.1|) is 
well defined only if the following condition is fulfilled 



which follows from Eqs. (|3.26|) and (|5.3|) . We remark that if the dynamics of the system 
evolves in one spatial dimension, Eq. (jB.9|) is trivially verified and the transformation in 
(15. 1J) can in all cases be accomplished. 
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